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Abstract 

We compute by means of the Bethe Ansatz the boundary S matrix for the open 
anisotropic spin-^ chain with diagonal boundary magnetic fields in the noncritical 
regime (A > 1). Our result, which is formulated in terms of g-gamma functions, agrees 
with the vertex-operator result of Jimbo et al. 



1 Introduction 



The concept of boundary S matrix in 1 + 1 dimensional integrable quantum field theory 
was precisely formulated by Ghoshal and Zamolodchikov in Ref. |J. They also developed 
there a bootstrap approach for determining such S matrices. Boundary S matrices can also 
be computed for integrable quantum spin chains by a direct Bethe-Ansatz approach which 
was proposed in Ref. ||. [] This method has been used until now only for isotropic models 
01) " @- We recently simplified this method in Ref. [[|. In the present Letter, we take 
advantage of this simplification to analyze an anisotropic model, namely the open XXZ spin 
chain 



4 n=l 

+ sinh 7] coth(?7^_) a{ + sinh r\ coth(r)£ + ) o z N \ , (1) 



where the real parameters £± > 1/2 correspond to boundary magnetic fields. For simplicity, 
we restrict our attention to the case A = cosh 77 > 1, which corresponds to the noncritical 
regime in which there is a nonzero gap. (See e.g. Ref. |10[|.) Our result for the boundary 
S matrix, which is formulated in terms of (/-gamma functions [TTJ, agrees with the result 
found by Jimbo et al. by means of the vertex operator approach. In the limit rj — > 0, 
we recover the results of Refs. 0] and H] . 



2 Bethe Ansatz and one-hole state 



In this Section we review the exact Bethe Ansatz (BA) solution of the open XXZ chain, and 
we compute the root and hole density for the BA state consisting of a single hole. 
The eigenvalues of Tt and S z = lJ2n=i a n are given by |L3 
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E = — sinh 2 77 7 r\ 7 ^ , (2) 

2 a= i sin r] (A Q - I J sinr/ (\ a + |J 

S z = -~M, (3) 



1 This method is a generalization of the approach developed by Korepin-Andrei-Destri ||, Q to calculate 
bulk two-particle S matrices. See also M. 
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where Ai , • • • , Am satisfy the BA equations 

sin 77 \X a + - -)) sin 77 (A a + - -)) / sin 77 (\ a + §) 



2JV 



sinr/ (A Q - £(£+ - |)J sin 77 (A Q - z(£_ - |)J ysinr/ (A a - 



smr/(A a - A^ + 2)smr](A a + A / 3 + z) 

= 11 ~ 7\ \ T~ 77 — n w> a = l,---,M. (4) 

^ an 77 (A Q - A^ - 1) sin 77 (X a + A^ - z) 

(We neglect in Eq. (H) additional terms which are independent of {A a }.) 
Introducing the notation 

sm V (x + f) cosr/fA + f) 

e n (\) = '-) 2 -L , g n (\) = Ur , (5) 

1 ' sinry(A-f) 1 ' cosr? ( A -f) 1 ' 

the BA equations take the more compact form 

M 

e 2 5+-i(A Q ) e 2€ __i(A a ) gi{K) ei(A a ) 2JV+1 - - e 2 (A Q - Xp) e 2 (\ a + Xp) . (6) 

(3=1 

Note that the factor gi(X a ) is absent in the isotropic limit 77 — > 0. 

Without loss of generality, we restrict 77 > 0. Moreover, the requirement that BA solutions 
correspond to independent BA states leads to the restriction (see || and references therein) 
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Following |§, we focus now on the BA state consisting of a single hole. This state lies in 
the sector N = odd with M = y — ~ and {A a } real. This state has S z = +|. 

Since Eq. (||) involves only products of phases, it is useful to take the logarithm. In this 
way we arrive at the desired form of the BA equations 



'a ; 



h(X a ) = J 

where the so-called counting function h(X) is given by 

h(X) = ^{(2N+l)q 1 (X)+r 1 (X) + q 2 ^ 1 (X)+q 2 ^ 1 (X) 

M 

Efe(A-A /3 )+g 2 (A + A /3 )]}, (9) 



p=i 
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g n (A) and r n (A) are odd monotonic-increasing functions defined by 

-tt < q n (X) < n 



?n(A) = tt + zloge n (A) 
r n (A) = ilogp n (A), 



-TT 



< r n (A) < 7r 



(10) 
(11) 



and {</«} are certain integers which serve as "quantum numbers" that parametrize the BA 
state. (See e.g. p6|| .) 

We shall need in the next section the root and hole density cr(A) for the one- hole BA 
state, which is defined by 



a(X) 



1 dh(X) 
N dX 



(12) 



To calculate this quantity, we must pass with care from the sum in h(X) to an integral. 
Indeed, with the help of the Euler-Maclaurin formula for approximating sums by integrals, 
and using the fact that the solutions A = , ^ of the BA equations are excluded, one can 
derive?] for a state with v holes the following general result: 



1 ^ r— 1 v 1 T 7r 

jf X>( A «) = / 2 " dX a W - w - ^ 0(0) + g(^) 



(13) 



(plus terms that are of higher order in 1/N), where g(X) is an arbitrary function, and {X a } 
are the hole rapidities. 

Using the above result, we obtain the integral equation 



a 



(A) = 2ai(A) + j^{aiW + h{X) + a 2 (A) + 63(A) + a 2€+ _x(A) + a 2? __i(A) 



+a 2 (X~X) + a 2 (X + X)}~ f 2v dX'[a 2 (X-X') + a 2 (X + X')]a(X'), A>0 (14) 



where A is the hole rapidity, and 

1 d 



. „ r] sinh^n) 

a n{X) = ^-TT?n(A) = r~— j—r 

2tt aX 7r cosh(T]Ti) — cos{2r]X) 

l_d_ (X ,_V sinh^n) 
n[ } 2ndX Tn[ } TrcoshM + cos^A) 



a n (A ± 



2rj J 



2 The argument is a slight modification of the one presented in the appendix of Ref. 



(15) 
(16) 
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Note that the b n (X) terms are absent from the integral equation in the isotropic limit. The 
b\ term originates from the factor g\ in the BA Eqs. and the 62 term originates from 
the last term in Eq. QIB|). 

The symmetric density o s {X) defined by 

N f a(X) X > 
°^ = {e(-\) A<0 < 17 > 

can now be readily found with the help of Fourier transforms, for which we use the following 
conventions: 

f(X) = - E e- 2i " fe V(A0 , /(*) = / dA e 2 ^ A /(A) . (18) 



k=—oo 

Indeed, we find 



a.(A) = 2*(A) + ir«(A), (19) 

where 

r (+) (A) = s(A) + (A) + J( A) + L(X) + (A) + J M (A) + J( A - A) + J( A + A) , (20) 
and 

5 = -^, J=^V, ji +) = ^, K = -^, L=-^, (21) 
1 + a 2 1 + a 2 " 1 + a 2 1 + «2 1 + «2 

with 

a n (fc) = e-^' fc l , b n (k) = (-) k a n {k) . (22) 
Note that the Fourier series for J±~\X) converges for £± > 1/2. 

3 Boundary S Matrix 

The boundary S matrix has the diagonal form 

°i ). (23) 
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The matrix elements a(X , £±) and (3{\ , £±) are the boundary scattering amplitudes for one- 
hole states with S z = +| and S 12 = — ~, respectively. 
We first compute a(A,£±). Setting 



a(A,£ ± )= e ^±\ 

we obtain by a calculation completely analogous to the one in Ref. || the result 
$«(A) = 0(A,£+) + 0(A,£_) = 2tt / A r (+) (A) dA + const. 

JO 

Recalling the result ( PU[ ) for r( + )(A), and using the fact 



(24) 



(25) 



J (A - A) + J(A + A) 



dX 



2J(2A) dA 



we obtain 



0(A , C±) = / s(A) + K(A) + J(A) + L(A) + 2 J(2A) + 2 4 +) (A) 



(26) 



(27) 



We now use Eqs. fllq) , (pl|), (p2|) to explicitly write the integrand as a Fourier series. 
Performing the A integration, using the identity 



oo ^—2r)kx 



J 1 J- 
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log 
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r> (^) 



i 



log(l - q 4 ) 



(28) 



where q = e v and r ? (a;) is the g-analogue of the Euler gamma function (see Appendix), 
and using also the g-analogue of the duplication formula [11] 



r q (2x) r>(~) = (1 + q) 2 *- 1 i>(x) + ~) , 
we obtain the following result for a(A,£±) (up to a rapidity-independent phase factor): 



(29) 



a(A,£i 



-4iA 
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r> (^ + i( 2 e ± -iy 


) r> (f + ±(2£± + 1)) 


t q 8 (t + i) 


r> + 1 


) r>(f + J(2e±-i)) 


r> (# + 1(2^ + 1)) 



.(30) 



We turn now to the computation of /3(A , £±), for which we must consider a one-hole state 
with S z = — |. Instead of taking the pseudovacuum to be the ferromagnetic state with all 
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spins up as we have done so far, we now take the pseudovacuum to be the ferromagnetic 
state with all spins down. The expression (fj) for the energy eigenvalues remains the same, 
the expression (|3[) for the S z eigenvalues becomes 

S* = M-j, (31) 



and |To] there is a change £± — > — £± in the BA Eqs. (|J). 

We focus on the BA state consisting of one hole (M = — — ~ with {\ a } real), which 
evidently has S z = — |. The corresponding function n -*(A) is the same as r^ + \\) (see Eq. 

(|20|)), except that now J±~\\) is replaced by J± \\), with Fourier transform 

J± -~TT^- (32) 



We observe that 



ff(A,£-) ff(A,£+) _ e i 2 7r/ o X [r(-)(A)-r(+)(A)] dA 

a(A,£_) a(A,f+) 



(33) 



Using the identity 



we conclude that 



4" , (A)-4 +, (A) = -a 2 £ ± -i(A), (34) 



%T\--^W- (35) 

We have fixed the sign in Eq. ( |3~5|) by demanding that /C(A , £±) be proportional to the unit 
matrix for A = 0. 



4 Discussion 

Our final result for the boundary S matrix of the XXZ chain is 

K(U ± ) = a(A,« ± )(J -eJ_ l( A))< < 36 > 
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where a (A , £±) is given by Eq. (|30|) . It can be shown that this result agrees with the one 
found by Jimbo et al. |12| by means of the vertex operator approach. In the isotropic limit 
rj — > 0, we see that q — > 1 and T q (x) — > T(x); hence, we recover the boundary S matrix of 
the XXX chain 0, §. 

It would be interesting to see if this Bethe Ansatz method can be extended to the critical 
regime \q\ = 1, which is outside of the domain of the vertex operator approach. 

The R matrix for the XXZ chain is associated with the fundamental representation of 
The present work opens the way to calculating boundary S matrices for spin chains whose 
R matrices are associated with the fundamental representation of any affine Lie algebra. For 
higher representations, the ground state involves complex strings, and the analysis is more 
complicated. We hope to address these and related questions in the near future. 
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6 Appendix 



Here we prove the identity 



fc=i - 1 



e ~2r)kx j 

+ e" 2 '^ k 



log 



r> (f ) 
IV (¥' 



log(l-g 4 ) 



(37) 



where q = e v and T q (x) is the g-analogue of the Euler gamma function, which is defined 



11 as 



r 9 (x) = (i-g) 1 ^n 

j=0 



It is convenient to first consider the sum 



< q < 1 



(3* 



oo g— 2r/kx 



S(X) ~ £ 1 + e- 2 ^ ' 



(39) 
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Expanding the denominator in an infinite series and then interchanging the order of sum- 
mations, we obtain 



where 



S(x) 



n^—2r\kn 



^ e -2r]kx Y 

k=l n=0 

oo oo 

^^(_) n ^-2r\k(x+n) 

n=0 k=l 

-2r){x+2m) 



-2rj{x+2m+l) 



| 1 — e -2rj(x+2m) ^ _ e -2r){x+2m+l) 

1 



log q 4 



^ 4 (I) ~ ^ ( 



x + 1 



^iog r g (x) 



oo n x+j 



log(l -g) +\ogqJ2 — 



x+j 



(40) 
(41) 

(42) 
(43) 

(44) 
(45) 



Integrating the result fl4"3j) with respect to x, and evaluating the integration constant by 
considering the limit x — > oo, we obtain the desired identity Eq. (]37f) . 
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